UNCLASSIFIED 


i 


fiepAoduced 
k4f  tJte 


ARMED  SERVICES  TECHNICAL  INFORMAnON  AGENCY 
ARUNCTON  HALL  STATION 
ARLINGTON  12.  VIRGINU 


UNCLASSIFIED 


IfOTICK:  When  government  or  other  drawings,  speci¬ 
fications  or  ether  data  are  used  for  any  puzpc:!e 
other  chan  In  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  idiat soever;  and  the  fact  that  the  Govern¬ 
ment  nay  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  la  not  to  be  regarded  by  Implication  or  other¬ 
wise  as  In  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rl^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  Inventlcn  that  nay  In  any  way  be  related 
thereto. 


CATALOGED  BY  ASTIA 


8  CO 


Contract  No  af  49,638  -30^ 
Projict  No  9762 


BlOrS  VARIATIONAL  PRINCIPLE 
IN 

HEAT  CONDUCTION 


Thonvu  J.  |jirdn«r  «nd  Fr«d*rki  V.  PoN* 


lOLYTECHNU.  INSIITUTt  OF  BROOKLYN 


OEPARTMfNT 

«• 

AERt^SPACE  INCilNEERINi; 

Anti 

APPIIEI)  MU  MANU'n 
JULY 


PlUL  REPORT  NO.  59;^ 


AFOSR  1275 


BIOT'S  VARIATIONAL  PRINCIPLE 
IN 

HEAT  CONDUCTION 


bv 


Thomas  J.  Lardner  and  Fradarlck  V*  Pohla 


Polytechnic  Institute 
Department 
Aerospace  Engineering  and 


of  Brooklyn 
of 

Applied  Mechanica 


July  1961 


Pibal  Ranert  No.  B95 

This  research  was  supported  by  the  Air  Force  Office  of  Scientific 
Research  of  the  Air  Research  and  Development  Command  under  Ccn<^ 
tract  AF  *»o(63S).30i}. 


li 


bibliographical  coiitrol  sheet 


1: - OrlqlnatlPQ  ao*ncv  and  mcnitorlno  «Q»ncv« 

O.A.i  Polytechnic  Institute  of  Brooklyn,  New  York 

M.A.t  Air  Research  and  Development  Command,  Air  Force 
Office  of  Scientific  Research 

Zi _ Orlqlnatlno  aoencv  and/or  monitoring  aoenc/  report  numbers 

O.A. t  PIBAL  Report  No.  595 

M.A.t  AFOSR  12/5 

It _ Title  and  classification  of  titlet 

BICT'S  VARIATIONAL  PRINCIPLE  TN  HEAT  CONDUCflON.  (UNCLASSIFIED) 

lit _ Personal  authort 

Thomas  J.  Lardner  and  Frederick  V.  Pohle 


_ 

_ Date  of  reoortt  July  1961 

6, _ 

— iJLlll*  59 

7, 

Illustrative  material  1  2 

figures 

9i 

Precared  for  Contract  No.  i 

AF  99(638) -302 

' 

Preoared  for  Proiect  No. i 

9782 

10. 

Security  classification! 

UNCLASSIFIED 

lit. 

■  ■  Pit  \  1 1  But  u  blui;  up  t » 

In  accordance  with  enclosed  dls 

^rlbution  list. 

12t _ The  literature  cn  Biuv's  variational  principle 

for  heat  conduction  and  the  thermodynamic  foundations  of  the  prio¬ 
ri  pie  are  reviewed.  An  additional  example  to  those  presented  in 
a  prevlou*  paper  Is  given.  Tnls  example  treats  the  heating  of  slabs 
exposed  to  time-dependent  heat  fluxes  and  a  specific  example  of  a 
triangular  heat  pulse  Is  presented  In  detail* 


TABLE  OF  CONTENTS 


Page 

ABSTRACT  iv 

SECTION  Ii  LIST  OF  SYMBOLS  1 

SECTION  III  INTRODUCTION  4 

SECTION  nil  LIIERATURF  REVIEW  7 

SECTION  IVi  THERMODYNAMIC  FOUNDATIONS  OF  THE  VARI-  18 

ATIONAL  PRINCIPLE 

SECTION  Vi  THE  HEATING  OF  SLABS  EXPOSED  TO  TIME-  3"^ 

DEPENDENT  HEAT  FLUXES 

SECTION  VIi  CONCLUDING  REMARKS  50 

SECTION  Vllt  REFERENCES  51 

FIGURES  ^  56 


ABSTRAHT 


The  litereture  on  Biot's  variational  principle  for 
heat  conduction  and  the  thermodynamic  foundations  of  the  prin¬ 
ciple  are  reviewed.  An  additional  example  to  those  presented 
in  a  previous  paper  is  given.  This  example  treats  the  heating 
of  slabs  exposed  to  time-dependent  heat  fluxes  and  a  specific 
example  of  a  triangular  heat  pulse  is  presented  in  detail* 


SECTION  I 


LIST  OF  SYMBOLS 


constant,  in  entropy  expression  for  neign borino  tner"i''“ 
dynamic  state,  Eq.  (  1  ' 

slab  thickness 

thermodynamic  coefficients  in  relations  between  forces 
and  fluxes,  Eo.  (4) 

spe  c  f :  c  ne  ac 

volume  element 

0  i  s  s  :  pa  ;  on  function 

operation ai  form  o:  dissipation  function 

V. 

exponential  function  of  x 
defined  by  £q .  ^  26 / 
defined  by  Eq .  ' 25 ) 

neat  transfer  coefficient;  heat  flow  fiom  System  I  to 
System  !  I 

heat  flow  from  energy  reservoir  m  System  I 
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surface  temperature 
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generalized  coordinate 
thermal  force 
generalized  thermal  force 
entropy:  entropy  of  System  1 

ij 

entropy  of  total  system 

defined  by  Eqs.  (131 

defined  by  Eqs.  '13) 

entropy  in* low  to  System  I 

entropy  of  constant  temoerature  reservoir 

entropy  of  energy  reservoir 

time 

pe-'etriticn  time 
temper  a  t  u  r  e 
reser.'Cir  temperature 
interna’  energy 
ae f  : ned  by  Eqs .  ^  ■ 3 ; 
ae ♦  .  ned  cy  Eqs .  '  i 3 ) 

♦nermai  potential  function 

d  e  f  1  n  o ,  j  i;  y  Eqs.  ;  i  3 ) 

defined  ov  Eqs.  13) 

spatial  coordinate 

generalized  thormodynam i c  force 

defined  by  £q.  (26) 

tnerr-’l  diffusivityt  k/c 

variation 
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dimensionless  penetration  depth 

temperature  difference  measured  from  oqu  i  1  Ibriufr  ’ 
temperatuie 

dimensionless  time 

•j 

aimensionless  penetration  time 
dimensionless  surface  temperature 
surface  temperature  parameter 
d .  dt 

i  (  '.3 X  )  1  (  3 / 9 V ')  e  n  <07  ) 
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SECTION  II.  INTRODUCTION 


The  high  speeds  of  niodeTn  aircraft  end  guided  mis- 

j 

siles  have  introduced  new  proclems  into  the  design  of  struc¬ 
tural  elements  for  such  vehicles.  These  problems  are  a  conse¬ 
quence  of  the  large  heat  input  to  the  structure  arising  due  to 
the  increase  in  temperature  of  the  air  in  immediate  contact 
with  the  surface  cf  the  veniule.  Th.s  increase  in  temperature 
is  caused  either  bv  shocK  compression  at  the  blunt  leading  edi^a 
or  by  friction  In  tpe  boundary  layer.  Similar  type  of  heat  in¬ 
puts  alec  crise  c^e  to  the  flow  of  nigh  energy  gases  in  rocket 
motors  and  nucleai  reactors  although  the  source  of  heat  is 
ohvsicallv  cifferent.  In  all  cf  these  situations  the  tempera¬ 
ture  ri'e  in  the  stiucturai  components  must  oe  accounted  for 
ir  the  ac-sign  of  tneentire  structure.  This  requires  among  otner 
prociems  the  de te rm i n at i cn  from  a  neat  *onduction  analysis  of 
tne  temperature  aistrinuticn  in  tne  structure  under  the  known 
thermal  ’ npu  t  s . 

Howeve-,  tne  large  neat  inputs  and  the  accompanying 
high  temperatures  will  necessitate  a  consideration  of  the  tem- 
peraluie  dependence  of  the  thermal  properties  of  the  heated 
oody.  The  introduction  of  these  two  fffer.ts  will  make  the 
gov er fling  neat  conduction  equation  ^nd  noundary  conditions  non¬ 
linear,  a  problem  which  can  rot  be  treated  by  the  usual  classi¬ 
cal  methods  of  neat  conduction  analysis.  Refs.  [l.  12]. 


In  such 


s^tu^itions,  3ppr3>imate  aralvtic^l  r.  r  numerical  metnods  ^:p 
required  in  ordpr  to  obtain  «  lolution.  A  discu&v.io''  cf  somo 
:f  tbp  ^nalvt:'-2l  and  numer;'*’.  iretb.ods  of  solution  was  pre 

epntod  in  Reis.  ['.!<*  i”'*-  ‘^J*  ■' 

One  :♦  tr.p  most  important  cf  tbP  ana  vtical  metnods 

fc-  npa‘  -  c''n;>  - 1 cr.  is  tbe  variational  principle 
:-.ri.d..-ci  fcv  B.ot.  =^Pf<.  [2  to  ..]•  B'-ct  r.iic  iipproacbed  the 
crr.r.-rir  ♦rerr.  a  f  .nd..ient..  nt^ndc.  .nt  is  a  •ariaticnal  principiP 
t.-.prrc-^  .ist.-itv  wr.iir.  treat-  tne  lOL.p.ed  elastic  and  tner- 

»r-^****j  iCiriCT^Cl* 

r  a  ,  s3*r^:t-.va  *i  j 

-.ar-.a*  lOna;  tn^oren^  ct  -.ast-.citv  are  onvained  from 
tnP  .ar.at::na.  :  orn- 1  .  t .  cn .  If  tn^  Plast.t  effects  are  ig- 
:r*i.  -,^r:at.cr^:  -errr  c‘  -ne  neat  conauction  equation  is 

t  r  .  j  .  p  .1 . 

”  ^  s  pr.^ric.®  wa*  re.  iPwpn  ^rerr  a  tr  at  nemat  i  c  a  .  view 
*  cr»a  •  Ca'-e  0*  ’’  p  A  *  ''ono-ttion  in  Ref. 

.  . .  .  ..o  .  n.c'ussPd  fer  tne  case  m 


fp  n'n  n.S'uSSPCl 

'  p  r  r  I  d  c  "  d  .  t  i  c  n  i  -  .  c- .  d  t  r  r  .  p  a 


*  *  i  u  1 .  A  rr.  p  t  n  0  a 


*-•  ,''.-ii''g  '.-.tn  Doundarv  cc''o.t;o^:  'Anitr.  r.  ad  •'c  .  been  v.re-.e 

DPfere  bv  tne  -ariatioral  pri-.itlp  w.,  .ntrodi^ced  and  the  vari- 
atic-ai  principle  was  then  app-ied  to  a  number  of  different  one- 
dimensionai  heat  conduction  prorlem..  The  prcc.emt  treated  in¬ 
cluded  nor-llneav  heat  f;u>.  bcundary  conaitions.  t  ime -dependent 
aerodynamic  neatlngt  and  the  heating  cf  bodies  witn  tpmperat^-.e 
dependent  material  prcperties.  The  results  Obtained  for  tnese 


6. 


examples  demonstrated  the  appl' i  c  aiii  1 1  ty  cf  the  variational 
princip.e  tc  prcbiems  with  specified  ne«i  fiux  boundary  con¬ 
ditions. 

In  the  present  paper,  whictv  comnlements  Ref.  [13]. 
the  thermodvnatr  *.  c  f  cui'.dat  ions  of  Biot*8  variational  principle 
are  reviewed  together  with  'previous  work  on  variational  methods 
for  heat  conduction.  Since  Biot's  ''arialional  principle  in¬ 
cludes  the  field  cf  thcrmo-elastic itv.  a  discussion  of  other 
viTiationa.  principles  obtained  from  Biot's  f  orrr.u.  ation  Is  also 
piesentod.  In  th.s  wav.  the  principle  as  applied  to  the  analy¬ 
sis  of  rest  r o-Hi  "t n*'  p-oblem*  :c  put  i"  the  proper  perspective 
as  -t  specia<  c-ste  cf  a  more  general  principle. 

Section  III  contains  a  review  of  t)»e  literature  on 
previous  werK  employing  variational  principles  for  heat  conduc- 
tic'  ;nii  on  Sicf  s  var  lationa,  principle.  Section  IV  reviews 
I'o  tnermcii vn^nj  1^  foundations  ct  tne  general  variational  prin- 
.c.®  a-'d  1'*  fc"*'  ter  tne  specia.  case  of  neat  rcnductlcn. 

Tne  meiiiuu  fu;  i»'viv.ding  neat  f.ux  ooundary  conditions 
that  wa«  intrcd-i-ea  in  Ret.  [13]  l*  applied  in  Sect'ion  V  to  the 
prcDl*-rr.  of  the  heating  of  a  finite  slab  with  an  arbitrary  time- 
depc ".'ier't  heat  input.  Tr.e  nesting  of  a  slab  bv  a  tixangular 
neat  ouise  iS'  the».  rtisrussed  in  detail. 


SPCTIGN  III.  LITERATURE  nZVIZ'.V 

A  discussion  of  variational  principles  in  thtimc- 
elasticity  and  heat  conduction  is  ofesented  in  this  section. 
Although  the  present  paper  it  '-.oncerned  mainly  with  heat  con¬ 
duction*  the  discuss  ion  of  the  variational  principles  in 
thermo-elasticity  presented  belcw  provides  a  basis  for  the 
discussion  of  Blot's  variational  principle*  Refs.  [2  to  il]. 
This  principle  is  applicable  to  the  general  field  of  thermo¬ 
elasticity  and  hence  includes  the  neat  conduction  analysis  as 
yai O’  Iht’  ccmpiete  formu.at;cn.  Tnerefoi*,  by  presenting 
the  discussion  it,  this  way*  the  variational  principle  for  heat 
conduction  alone  is  /lewed  in  its  proper  perspective  as  a  spe¬ 
cific  ca^e  of  a  .xoie  generai  prin'-ipie.  In  uddl^'ori  to  these 
discussions  referei.ce  will  be  made  to  other  approximate  methods 
of  solution  for  the  neat  conduction  equation. 

Variotlonol  Princi&.es  in  Tr.e  i  tr.o  *  E I  a  s  1 1 1.  i  t  v 

B.ot*  Refs.  [4.  has  introduced  a  varictional 

principle  which  has  oppiication  to  tr.c  field  '.f  thermo-elasti¬ 
city  and  in  fact  to  mor-  generai  tno’-modynamic  systems*  The 
field  of  thcrmo-e I  as t : : ^ ♦y  includes*  In  the  most  general  sense, 
the  effeutfi  of  cnupl.ng  between  the  t empe r a tu re  and  elastic 
fields.  It  is  this  ''anational  principle*  restricted  to  the 
soecific  field  of  heat  conduction,  that  will  be  discussed  in 


this  pap<‘r 


The  ge'eral  variational  orinciple  was  derived  by 
Biot  on  the  b<ibiii  of  t ne ;  ir, c  J  v  ar. :  c  a .  arguments  appliea  to  a 
physical  system.  The  state  of  tnis  svstem  was  defined  bv  gen¬ 
eralized  state  .ariabies  and  Onsager  s  ;eciprocal  relations 
were  applied  to  obtain  Lagrang.an  equations  for  the  time  his¬ 
tories  of  tne  generalized  coc  r  vi  ;  n  a  t  e  s  .  Rets.  [?.  3.  This 

method  of  deduc.ng  tne  equat-o^s  will  oe  discussf-d  in  Settlor. 

IV.  Gnte  tne  general  equaticis  were  estaolished*  a  varia¬ 
tional  princip.e  eq-'va.ent  to  tnese  equations  was  introduced. 
Inis  V  r  ;  a  t ;  Cud  1  pwncicie  wan  applied  the  general  case  of 
thermo-el ast i c itv  ment.oned  aoove.  Refs.  [4,  ana  to  the 

case  where  tne  elasti*  effects  were  net  considered.  Refs. 

[5.  6.  9].  The  latter  case,  of  course,  a pc  lie  to  the  field 
of  heat  conduc  t  i  0" .  An  additior.«ji  field  of  application  of  the 
general  ■•ariaticnai  principle  which  r.as  been  treated  by  Biot 
.s  V  i  s  c  ce  1  a  s  1 1  c  1  tv  .  Refs.  [  1 1. .  .  I"  what  follows  in  this 

“'•'ion  r.cwc'.cr.  the  pri'cip.*.  wi..  oe  discussed  am  iG’atlon 
to  t  he  rmc  •  e  i  a  s  1 1 ;  t  V  wni:h  i-'ludrs,  i *ne  general  cense,  tne 
field  c*  neat  conduction. 

Biot  nas  sncwi  tnat  tne  gc  hg  eonatiens  for  equi- 

lib.lum  of  a  thermo-e 1 ast 1 c  solid  can  be  expressed  in  terms  of 
a  variation'!,  pr.nciple  under  specified  rp''strair.l  conditions. 
Ref.  [4],  I"  p.'.  r  t  ;  c  u .  ar  .  jf  the  equatiun  of  state  together  with 
conservation  c‘  energy  is  im'ssed  on  the  . ariational  "unction. 


9, 


th®  equilibrium  equations  and  Fourier-s  law  of  conduction 
aro  obtained  as  the  Euler  U  i  f  f  er*»nt  i  a  i  equations.  Tr.e  latter 
equation  with  energy  conservation  yields  the  coupled  neat  con¬ 
duction  equation.  This  principle  analogous  to  the  principle 
of  Minimum  Potential  Energy  m  Elasticity.  Refs.  [15.  16]»  in 
which  the  stress-strain  relations  are  imposed  and  the  equilib- 
brium  equations  are  tne  corresponding  Euler  equations. 

A  Cornpl  erne- tary  Energy  Principle  was  formulated  by 
Herrmann.  Ref.  [17],  and  leads  to  the  stres s -displ acement  equa¬ 
tions  together  with  the  conservation  of  energy  expression  under 
the  asi -iiipL  1  v:;  ci.dt  ihe  varieo  stress  ana  temperature  fields 
satisfy  tne  equilibrium  equations  and  Fourier'-s  law  of  conduc- 
ihe  emphasis  on  Fourier's  law  of  conduction  is  not 
clearly  stated  by  Herrmann?  this  point*  however:  is  particularly 
important  sii-ce  it  forms  the  basis  of  the  application  of  the 
variational  principle  In  heat  conduction.  This  was  noted  In  Ref. 
[13].  It  IS  seen. from  tne  above  that  this  principle  is  analagous 
to  tne  principle  of  Minimi'm  Complementary  Fnergy,  Refs.  [l5.  16]« 
in  which  the  eq'ii  librium  equations  restrict  the  arbiv.rA^y  stress 
variations  to  yield  the  compatibility  equations. 

Herrmann,  in  two  later  papers.  Refs.  [18.  19].  has 
considered  an  extension  of  Helssner-s  variational  principle* 

Refs.  [20.  2*.  22].  into  thermo-elasticity  in  which  all  the  field 
equations  are  obtained  as  a  consequence  of  a  variational  principle. 
Tne  principle  has  been  outlined  fox  the  case  of  a  uniaxial  state 


10. 


of  stress*  Ref.  [18],  and  generalized  to  a  three-dimensional 
anisotropic  oody.  Ref.  [19]. 

However*  this  extension  to  a  more  general  varia¬ 
tional  principle  introduced  an  additional  variable  called  the 
"thermal  dis-equ’ 1 ibrium  force*  which  has  not  been  employed 
by  Biot  or  by  Herrmann  in  setting  up  the  separate  variational 
principles  for  thermo-e 1 ast i c itv  mentioned  above,  "his  force 
has  the  form  used  i*’'  irreversible  thermodynamics  since  it  is 
related  to  the  gradient  of  the  temperature  field.  Moreover, 
it  provides  the  connecting  link  between  the  law  of  conduction 
and  the  energy  eauation.  whicn  leads  to  the  coupled  heat  con¬ 
duction  equation. 

With  the  introduction  cf  tnis  independent  va’-’able, 
the  governing  field  equations  were  octaineo  from  the  function¬ 
al  expression  given  by  Herrmann  in  a  manner  similar  to  that 
outlined  bv  Heissner.  Comparison  with  neighboring  states 
which  do  not  satisfy  the  variaticnal  principle  showed  that  the 
principle  is  th.’t  of  a  static narv  va.ue  problem  rather  than. a 
true  maximum  or  ...inimum;  again  this  paralleled  the  development 
presented  bv  Meissner. 

An  additional  formulation  of  a  variational  principle 
to  include  thermo-elasticity,  plasticitv  and  creep  deformation 
has  been  presented  by  Besseling.  Refs.  [23,  .2^].  This  vari¬ 
ational  principle  was  stated  in  terms  o^  the  physical  displace¬ 
ment  and  entropy  displacement  fields  and  leads  in  the  general 


1  1 
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t:?se  to  a  stationary  value  p^'-nciple.  Tne  entropy  displacement 
w-js  tHp  same  as  that  used  bv  diet-  Ref.  [4],  Furthermore,  an 
appli.;atIor  cf  this  variational  principle  to  a  problem  of  struc- 
tural  dampirg  vvas  presented. 

Variation, a'  Principles  m  Heat  Conduction 

Weirer  Ret.  [25].  presented  a  metnod  of  soltition  for 
the  heat  cc'ducticn  equation  which  combined  the  standard  Laplaue 
tech'll  que  witn  that  cf  tne  Galerkin  technique  (see.  for  example. 
Refs.  [.5]  and  [i6j'.  In  this  ♦' ;-niu  1  at  i  on ,  the  transformed  heat 
-onducticn  rq-utich  in  the  Lapl-.e  was  solved  approxi- 

matelv  ov  tne  Galerxin  technique.  Tne  approximate  solution  ob¬ 
tained  m  this  manner  was  then  transformed  back  tn  recover  the 
time  variabie.  The  conclusion  drawn  by  We.ner  In  applying  this 
method  to  sene  simple  examples  and  to  the  flew  of  heat  in  an 
angle  seci.uu  was  mat  for  prcolems  cf  phvsital  interest  consider- 
>r. »  T'-p-  tat:C''3i  labor  i lequirrd. 

W-iir.  .lu.  Ref.  [26]  n-is  presented  a  variational  prin* 
iic.e  A,  •'ll!"  i'l  roiuced  a  tem.po'^itu’P  field  that  Wvas  the  mirror 
im^gp  z*  tne  pr.vsica.  temperature  d  i  a  t  r  i  but  ic"  .  This  artificial 
•.  emp  e  :  i  tu  :  e  .<,'■  he  i  n  ter  pr  e  t  erj  as  fit  corresponding  to  a  negative 
theimai  '  c”duc  1 1 V  1  ty  .  T  m^tiicd  was  app.ieo  ry  assuming  a  pciy- 
''cmial  temperature  distribution  to  oe  used  in  the  voriatvonal 
fu“cticri  ficm  which  the  coefficients  of  the  polynomial  were  found. 
An  example  in\civlrg  the  neating  of  a  finite  slab  with  fixed  sur- 


face  temperatures  showed  poor  j*greeir.er.t  with  the  ejiact  solution 
fcr  a’’'  appreciable  t  '.me  range,  m  particular,  for  short  times. 

Roser,  Ref.  [2~].  formulated  a  variational  principle 
based  upon  a  rr  od  i  f  1  :  a  1 1  on  of  Onsarer-s  principle  of  minli.,ti'n  dis¬ 
sipation.  The  functional  expression  introduced  was  expressed  In 
terms  of  tne  temperature  gradients  and  the  time  derivative  o- 
tne  temperatvre.  This  functional  expression  was  varied  with  res¬ 
pect  to  the  temperature  witn  the  time  derivative  held  fixed. 
Appropriate  boundary  conditions  were  eivfcrced  so  that  the  heat 

'  o-^du  1 1 1  c ''  equation  was  cotai'^ed  from  the  /ariatioral  fui.ction. 

■ 

innile  '.ne  .ariat'ora.  princip  e  is  'insistent  within  the  frame- 
wcr?:  prese-^ted  a  .ogira.  inconsistency  does  exist  If  conslder- 
I'.ic''  IS  gi.-e'-  to  c  c  n  se  r  /  a  •.  n  n  of  energy.  Mo  applications  of 
thic  rr:'r:p,P  we;e  presented  cv  Rosen  so  tnat  no  statements  as 
tc  app 1 1 • ab i 1  .  tv  rf  tne  prir.icle  : an  be  made. 

A  •jriationai  pT.'^cipie  aimost  opposite  to  that  presen¬ 
ted  cv  •<cs«=-n  .vds  gi  -en  pv  Chambers.  Ref.  [28]  In  this  varia- 
tio’“-iI  tJ  ’  •'  ;p:'  bds®d  upc'*  Ww.f.  oi  "erivel.  Ref.  [2^].  a  vena- 
tiona.  function  was  c*p-*‘se«»d  in  termc  of  tne  temperatur®.  the 
time  rate  r,f  hange  of  the  temce  r  xf; -e .  me  temperature  gradient 
and  the  heat  flux.  The  var.ation  was  performod  with  reapect  to 
the  time  derivative  of  the  temperature  whMe  the  temperature  It¬ 
self  remained  unchanged.  However,  although  the  heat  conduction 
eguatinn  and  tne  appropriate  boundary  conditions  were  obtained 
in  ^  fcrm?i  manner,  a  physical  inconsistency  exists  In  this  prlr.- 


ciple.  This  is  due  to  the  introduction  of  the  heat  flux  with¬ 
out  '•.onsidering  its  relation  to  the  time  rate  of  change  of  the 
temperature  through  tne  expression  for  energy  conservation. 

N’o  illustrations  of  the  principle  are  presented. 

Tne  above  two  variational  principles  illustrate  the 
.  '  or.  ?  i  i  t  :  i  es  that  ma,'  v.-.en  a  variational  principle  Is 

'.rticiuced  without  .*1*1  underlying  physical  concept.  In  parti- 
'ulir,  the  3e>,ond  principle  is  rot  correct  frem  a  ohysical 
cr.”t  of  .:ew.  On  the  other  hand,  the  variational  principle 
-;''trodu;ed  cv  Blot  is  ir.ctivated  by  strong  physical  reasoning. 

Tn :s  prin.ip.e  which  was  discussed  above  has.  in  addition*  the 
aavantage  ot  general  it/.  However,  as  was  mentioned  before* 
f'e  /-iriational  principle  will  be  applied  in  this  paper  to  heat 
corjjction  processes  cniy.  iins  is  equivalent  to  neglecting 
the  rlnsti-  effects  in  the  thermo -e 1  as t i c  variational  principle. 

The  mol  ivat  i  cn  foi  the  principle  from  thermodynamic 
f'^icer^t:  c'-s  will  be  presented  in  Sectio"  IV.  An  additional 
iiicz..3SiOh  cf  the  the.iTiOdVMorrtic  p**". Cip.  —  s  in  oiOt^^  Lneuiy  won 
fc-nj  in  Ref.  [30] 

The  vari atiorai  principle  was  discussed  from  a  mathe- 
maii.al  ooint  of  v.ew  ir.  Ref.  [13].  That  is.  the  mathematical 


steps  snowing  tne  equivalence  ot  tne  variational  principle  to 
th*  neat  conduction  equation  was  piesented.  Similar  discussions 
of  :.he  variational  principle  were  given  in  Refs.  [31,  32], 
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In  addition,  tht  problta  of  tho  dotomination  of  the  teaperature 
dial  ;.bution  in  a  flange-web  coabination  (angle  section)  which 
was  treated  by  Bioty  Ref.  C^3t  war  further  discussed  in  Refs. 

[33,34].,  An  extension  of  the  variational  aethod  for  this  pro- 
blea  to  two  diaensions  was  given  by  Levinson,  Ref.  C3&]*  These 
references,  with  the  exception  of  Refs.  [13, 3b],  do  not  present 
any  new  foraulations  or  results  of  the  variational  principle. 

Hnwever,  the  principle  was  aodified  for  the  ease  of  both  tran¬ 
sient  and  steady  forced  convection  heat  transfer  in  Refs.  [36 
to  39]. 

Reference  [37],  Part  I,  contains  a  foxaulation  of  a 
variational  principle  for  the  transient  heat  convection  equation 
analogous, to  that  foraulated  by  Biot  for  the  heat  conduction 
equation.  Gupta,  Ref.  [38],  has  extended  this  principle  for 
convection  to  the  case  where  the  medlua  is  anisotropic.  Both 
cf  these  principles  require  a  aodiflcatlon  of  Biot*s  work  to 
account  for  tne  motion  of  the  medlua  and  for  the  viscous  dissipation., 
If  the  latter  effect  is  neglected,  the  resulting  heat  convection 
equation  is  identical  to  the  heat  conduction  equation  for  a  moving 
medium.  In  this  case  the  dissipation  function  in  the  principle 
is  aodified  by  expressing  the  flux  condition  in  terms  of  the  velo¬ 
city  of  the  aedlumi 
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^9*'***^»  C373j  Part  II)  introducad  a  variational 

principle  for  staidy  haat  convection  in  channel  flow  under  the 
assuaption  that  the  effect  of  axial  conduction  can  be  neglected. 
Again  thii  variational  principle  is  siailar  to  that  of  Biot's. 
Agrawal,  in  a  later  paper)  Ref.  [36])  extended  his  previous  prin 
ciple  to  account  tor  the  vanishing  of  the  velocity  at  the  walls 
of  the  channel 

The  restriction  of  negligible  axial  heat  conduction 
ha,  been  set  aside  in  the  work  of  GuptS)  Ref.  [3»])  and  a 
variational  principle  has  been  developed  for  the  coeiplete  steady 
heat  convection  eoualien  for  channel  flows.  In  this  work)  as 
well  as  in  the  above  references  for  convection,  Lagrangian  type 
of  equations  have  been  forwulated  for  the  theraal  flow  field 
using  the  concepts  of  theraal  potential)  dissipation  function 
end  generallaed  force  The  use  of  these  concepts  clearly  shows 
the  rwl4ti.on  of  this  work  on  convection  to  the  work  of  Biot  In 
heat  conduction. 

Citron,  Refs  [eo,  *l],  has  discussed  Blot's  principle 
with  reference  to  ablation  problems.  In  CitronV  work  the  verie> 
tionei  principle  was  expressed  in  teras  of  e  functional  expres- 
•  lut.  the  heat  flux  end  the  teapereture  gradient.  The 

vertetion  was  carried  out  with  respect  to  the  heat  flux  with  the 
te»pt ratur#  gradient  held  fiwwH  This  particular  aethod,  however 


is  no;t  consistent  with  conservation  of  energy.  That  is,  any 
variations  in  the  heat  flux  aust  be  related  to  variations  in 
the  temperature  gradients  as  a  consequence  of  conservation  of 
er. . qy. 

Citron,  a»f.  [40],  also  presented  an  application  of 
Galerkin’s  aethod  (see  above)  to  the  ablation  problem.  This 
eethod  was  also  presented  in  Ref.  [32],  and  it  is  an  extension 
of  the  Galerkin  technique  to  the  case  where  the  approximate 
solution  1$  presented  as  any  function  of  the  time-dependant 
orbitrary  parameters  An  illustration  of  this  method  is  given 
m  Ref.  [32]. 

Other  Approximate  Methods  in  Heat  Conduction  Analysis 

The  rsmaiiung  portion  of  this  section  will  discuss 
additional  approximate  aathods  that  have  been  used  in  heat 
conduction  analysis  Although  numerical  methods  havo  been 
most  useful  for  the  determination  of  temperature  distributions, 
they  Will  no*  be  discussed  here.  However,  specific  numerical 
methods  were  discussed  in  Rof  [l3]  together  with  applications 
of  Biot's  variational  principla  to  prescribod  physical  pro- 
blams.  Hhile  these  numerical  methods  may  have  certain  advantage 
if  computational  assistance  is  eveileble,  they  roquiro 
in  any  subsequent  thermal  stress  calculation,  for  exampla, 
the  continued  application  of  numerical  methods.  This,  together 
with  the  temperature  problem  may  become  e  formidable  under¬ 
taking  In  addition,  if  an  analytical  solution  is  possiote, 


1? 


it  May  not  ba  in  a  convaniant  forn  for  nutterlcal  evaluation. 
TharafoT'ai  it  may  ba  desirable  in  these  casesj  for  the  purposes  of 
practical  coaputaticfti  to  obtain  a  siapla  approxiaate  analytical 
expression  for  the  teaperatura  distribution. 

Goodaan  has  introduced  an  approxlaate  aethod  called  the 
heat  balance  integral  technique  for  the  solution  of  the  heat  con¬ 
duction  equation!  this  aethod  has  been  u<ed  in  aany  studies  of 
heat  conduction  probleas,  (see  Ref>  [13]  for  a  list  of  references.) 
The  essential  idea  of  this  technique  is  to  satisfy  the  integrated 
heat  conduction  equation  with  an  assuaed  spatial  teaperatura 
piofile  that  aust  satisfy  the  prescribed  boundary  conditions. 

Many  of  the  results  obtained  in  Ref.  Cl3]  using  Biot's  varie* 
tional  principle  were  eoapared  «ith  this  heat  balance  technique. 

Green,  Ref.  has  presented  an  expansion  aethod 

for  a  general  fora  cf  the  heat  conduction  equation  which  is 
siailar  to  the  Galerkin  aethod^  An  application  of  a  aodifi* 
cation  of  this  aethod,  Ref  C32},  was  aentioned  previously  with 
reference  to  the  Gale-^'in  technique. 
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SECTION  1V=  THERMODYNAMIC  FOUNDATIONS  OF  THE 
VARIATIONAL  PRINCIPLE 

Th«  th«raodynaaic  foundations,  R»f»  [7],  which  lead 
to  the  general  fora  of  variational  principle  are  re¬ 

viewed  in  this  section<,  Aftei  Ihe  yutieral  fora  of  the  piin- 
ciple  is  developed,  the  discussion  will  be  restricted  to  the 
foundation  of  the  heat  conduction  variational  principle.. 

The  aajor  points  in  the  deveibpaent  of  the  general 
variational  principle  are  the  introduction  of  Onsager's  reci¬ 
procal  relations.  Refs.  (43  to  45j,  and  the  deteramation  of 
the  expressions  for  end  the  meanings  of  the  dissipation  func¬ 
tion  and  the  theraal  potential  Of  particular  laportanee  is  the 
^•'•fMi^ction  of  the  entropy  of  the  Total  syitea  in  toras  of 
the  entropy  and  the  internal  energy  of  the  primary  sub-syatea 
and  the  teaperature  of  the  heat  reservoir  Thu  leads  to  an 
expression  for  the  generelitetion  of  the  Helaholtr  free  energy 
for  s  syttea  at  a  nen-unifoia  teaperaturei  this  expression  is 
eslied  the  'Oenereitred  Free  Energy*  or  *Therasl  Potentiel** 
but  the  latter  tera  will  be  lised  ir  the  subsequent  discus- 
•  tens  T**e  equations  for  the  »tate  verieblee  of  the  ays  tea 
ere  first  developed  for  en  leoletod  systoa  end  then  for  a 
systoM  u-dor  extornal  toreos  Tho  typo  of  forcos  considorod 
mill  bo  of  a  thoraal  naturo  sinct  tho  finol  obltctive  is  to 
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apply  the  results  to  the  A^'^lysis  of  heat  conduction  probleos. 
However*  the  applicability  oi  the  wethod  for  the  deterwinatlon 
of  forces  ir  other  situatioiis  is  apparent  and  this  has  been  dis¬ 
cussed  by  Blot,  Refs.  [3,7]  Once  the  foraulation  for  the  systea 
has  been  ec  tablialieu  n:  sense  that  tne  Lagrangian  equations 

for  the  tiwe  histones  of  the  state  .'ariables  have  beer  found, 
a  variational  principle  equivalent  to  these  equations  it  ia- 
aediately  obtained  This  varia*io''8l  principle  restricted  to 
the  special  case  of  theres'.  f.elds  leads  to  the  heat  conduction 
equation  for  the  teaperatur#  this  wav  a  variational  princi¬ 

ple  aotivated  bv  tneraolvn#*::  arguments  can  p*  established  m  a 
consistent  aanrer  tor  neat  ionduetion  proces*e« 

Theraodvnaaic  Svsves 

The  entire  systea  to  o#  considered  consists  of  a  pn- 
aarv  subtyttea,  calUa  $».t#a  I,  to  which  a  heat  reservoir , 
called  Systea  11,  at  a  co'^start  teaperaturt  T^  ii  connected 
T**#  ert've  systew  vO^iisti'^  .j*  pptr  the  priaarv  systaa  and  the 
heat  reservoir  ;s  assuaed  to  be  isolated 

Injtialiv,  the  syite*  in  an  e^vilibriua  state  Ur, 
wich  the  values  e*  tne  internal  enerdv  are  taven 

as  equal  to  re^o  Deviatic**  froa  the  eqwt ubrsu*  stale  ere 
•vpressed  tn  teras  of  "  qene-iKrtd  coordinates  q^' which  ar* 
also  rero  in  i“'e  eqvilitriua  »ttt*  These  n  qenetalsred  co¬ 
ordinates  are  state  variables  oeftninq  the  state  of  the  e-tire 


system  at  any  instant  of  time.  The  valua  of  the  equilibrium 
temperature  is  and  the  local  deviat  ons  from  this  equili¬ 
brium  value  will  be  denoted  by  G;  these  deviations  are  con¬ 
sidered  as  state  variables  The  system  i  in  addition  assumed 
to  he  linear  in  the  sense  of  irreversible  thermodynamics 

This  requires  that  tne  flukes  in  the  system  are  linearly 
related  tc  the  corresponding  forces,  for  example,  that  the  heat 
flow  is  related  to  the  temperature  gradient  In  addition, 
production  of  entropv  m  the  system  must  be  written  as  a  pro¬ 
duct  of  tf  e  s'  d  the  corre  .ponding  forces  A  review 

or  these  bas.v  macros'-opic  .oncepts  tor  tne  thermodynamic 
theory  ot  irreversible  processes  was  given  bv  Miller  ir  Rtf 
f46].  Miller  discussed  the  f  a  t  tons  of  irreversible  thermo¬ 

dynamics,  wni.;n  ire  also  discussed  in  Refs  L^3  to  45l,  and 
prese'ted  the  basic  as^ .motion s  anj  motivation  of  the  theory 
’■he  system  at  'irsi  will  be  free  frcm  any  etternii 
applied  torces  Sin-e  tOr  e'^rrnpy  of  ^he  system  et  equi'i- 
hri’.m  ’'ds  a  maximum  vai'.'e  which  .i  t .  conven  i  en  t  ly  set  equa  I  to 
lern.  ihp  Bfitropy  c;  i  v  neighbvriny  state,  will  be  less,  than 
zero  I'his  value  c  entropy  can  be  e*V’...ressed  in^  terms  ot  the 
itate  '/itiable.  ’"hat  is,  the  entropv  at  jny  .  nei  ghfcbf  mg 

,tute  w  .,  th  •  •  i  te  var  i  »b  I  es  q  vahbt  w r  j  *  ten  i  n  ;  the  tor»‘ 
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»  c  • 
*  r' 


(1) 


where  S*  .s  the  entropy  ot  the  entire  sy«ite*,  the  4^^  are  the 
val-jes)  of  the  appropriate  derivatives  of  the  entropy  evalua- 
*teri  at  the  equiiibiiuffl  state  and  the  q^  are  the  state  variables 
The  summation  -onvent *0"  of  summing  over  repeated  indices  is 
Mseo  !»'  this  equating  and  in  the  equations  to  follow.. 

As  a  consequence  ot  the  deviation  from  the  equill- 
br  i>m  state,  r.«storing  forces  are  set  up  within  the  system. 

It  IS  assumed  that  these  forces  are  linearly  related  to  the 
time  rate  o*  change  of  the  state  variables,  the  fluxes,  q^ 

In  addition  i*  iS  as.sumed  that  the  linear  relations  satisfy 
tre  Onsagei  nypotheses,  Pefs  f ^3  to  4b]  This  requires  that 
tne  entiopv  product.on  m  rne  system  which  arises  as  the  sys¬ 
tem  -.s  returni-'g  to  the  equilibrium  state  be  written  in  the 
form 


T^fdS'/dt)  =  X^q^ 


(2) 


or  that 
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(3) 


This  latter  equation  states  that  the  forces  are 

equal  to  the  corresponding  gradients  of  the  entropy  in  the 

non -equiiibriua  state,  the  factor  is  introduced  for  con- 

ver.ierce  If  the  forces  X  are  now  e»oressed  in  terms  of 

1 

the  fluxes  q^  in  tne  following  form< 


y 


1 


Oneager’s  pr:!'::pl<‘  leads  to  the  result 


(4) 


Tnese  relations  intercorrect  the  various  thermodynamic  pro- 
:esses  taking  place  within  a  thermodynamic  system  and  they 
are  called  the  Onsager  reciprocal  relations,  R**fs.  [43  to 
4!^]  The  limitations  of  these  relations  are  also  discussed 
in  the  above  references  A  review  of  the  experimental  evidence 
tor  these  reciprocal  relations  lo  macroscopic  linear  relations 
IS  given  ;r.  Ref  [47]  The  experimental  evidence  given  indi* 
cates  that  the  Onsager  reciprocal  relations  are  satisfied  in 
cost  applications 
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In  view  of  Eq  ( 5) »  a  quadratic  function  D  exists  such 

that 

<1  =  ■),  - 

where 

D  =  ^ i/2)b^^  q^q^  (6) 

■^he  physical  significance  of  the  function  D  follows  iMediately 
from  Eq  (7)  which  e. presses  the  rate  of  entropy  production  as 
a  fv.  ret  ion  of  Llie  fuiuvs  d'‘d  flukes  q^» 

T^(dSVdt)  =  X,q^  =  q^q^  s  2D 

Therefore,  the  function  D  is  (l/2)T^  times  the  rate  of  entropy 
production  in  the  systewi  this  function  is  called  the  dissipation 
function  The  quadratic  form  of  D  expressed  by  Eq.  (6)  1»  then 
positive  oetinite  since  the  rate  cf  eutrepr  production  is  always 
positive 

In  addition  to  expressing  the  force  as  a  derivative 
of  the  dissipation  function,  the  entropy  of  the  system  Is  set 
eiqual  *o  a  function  V  called  tne  thermal  potential  such  that 


V 


(7) 


Therefore,  Eq.  (3)  together  with  Eq  (6)  can  be  written  in  the 
following  fornt 


(8) 


The<%e  equations  are  in  Lagrangian  forai  for  the  n  state  variables 

Tr.ey  are  analogous  to  the  Lagrangiar  equations  for  a  nechanlcal 
system  :n  which  V  is  the  potential  energy,  D  is  the  dissipation 
function  and  the  are  the  generalized  displacements.  The  ex-, 
t.ern^l  forces  acting  on.  the  mechanical  system  in  this  case  are 
equal  to  zero  however,  the  right  hand  side  o*  Eq.,  (3)  is  not 
zero  f  evterral  fcrr.es  aie  applied  A  similar  result  would 
be  e»pected  bv  analogy  for  the  thermal  system  under  external 
forces  However  before  starting  on  that  case  a  discussion  of 
thi»  t  ign  1  ♦  1  c  ar  ce  of  the  function  V  will  be  presented- 

Tne  Significance  of  the  thermal  potential  function 
V  ran  be  found  by  wiitiiig  the  entropy  of  the  entire  thermo- 
dynamic  system  ir  terms  of  the  entropies  of  the  two  subsystems* 

That  IS,  ♦he  entropy  S*  of  the  total  system  is 


whtrt  S  IS  th«  entropy  of  the  priwary  sub-system  and  is  the 
entropy  ot  the  constant  temperatute  reservoir.  During  the  pro- 

j 

cess  of  reaching  the  «quilibriua>  state^  heat  will  flow  in  general 
between  the  primary  system  and  the  reservoir.  If  this  heat  flow 
from  System  1  to  System  II  is  h,  then  conservation  of  energy 
for  System  I  gives 


U  =  -  h 


(10) 


where  U  is  the  internal  energy  of  System  T.  At  equilibrium 
both  h  and  U  are  eq’jal  to  zero  If  external  forces  are  acting 
cr.  System  I,  Eq  (10)  will  be  modified  by  additional  work  terms; 
this  will  be  treated  subsequently. 

Since  the  heat  flow  into  the  constant  temperature 
reservoir  is  Vnowr^  the  entropy  Sjj  is  found  in  terms  of  the 
intor^aX  energy  o*  the  primary  systemy  that  is. 


(11) 


Therefore,  from  Eqs  (^)»  (^)  and  (ll),  it  follows  that  the 
thermal  potential  function  V  can  be  written  in  the  form 
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Equation  (12)  is  similar  to  the  form  of  the  Helmholtz  free  energy 
for  the  System  1  except  for  the  oresence  of  the  reservoir  tempera¬ 
ture  Since  the  temperature  of  the  primary  system  is  not  speci¬ 
fied,  and  may  have  any  arbitrary  distribution,  this  form  has  greater 

applicability  than  the  usual  free  energy.  Blot,  Ref.  [2,7j,  calls 
this  function  the  grneralized  free  energy.  This  form  of  the  thermal 
potential  mill  be  used  later  in  the  discussion  of  external  forces. 

An  additional  form  of  the  f..,-ciion  V  in  lerms  of  tem¬ 

perature  can  also  be  found  and  it  has  a  more  convenient  form  for 
physical  application*;  than  that  given  above.  In  this  expression 
the  thermal  potenfai  is  considered  as  the  sum  of  two  partsj 
(U  a  value  obtained  when  all  the  state  variables  are  varied 
except  the  temperature  winch  is  held  fixed  anu  (2)  a  value  ob¬ 
tained  when  the  temperature  is  varied  and  the  remaining  state 
variaties  are  held  five-^  This  splitting  of  the  fvmction  V 
int-  two  parti  can  be  writtan  in  the  form 


V  = 


V  * 
r 


V 

c 


(13) 


\I  z  u-  r  s 

«  r  I  r 


y  z  u  -  T  s 
c  c  r  c 


where 
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These  latter  equations  e?«press  each  component  of  V  in  terms 
ot  th«  corresponding  values  of  the  internal  energies  and 
entropies  for  isothermal  changes*  and  and  for  con¬ 

stant  state  variable  changes*  U  and  S  That  is,  V  is 

O  *•*  e 

t.'ie  value  of  the  thermal  potential  for  an  isothermal  pro¬ 
cess  at  temperature  while  all  the  other  state  variables 
are  varied*  an  analogous  statement  applies  for 

The  term  can  be  evplicitlv  expressed  in  terms 

of  the  temperature  ♦  0  o*  System  I  This  value  of 

sor  the  entire  System  I  is 


(14) 


where  c  ’s  ffe  df  the  *if>*f.**ii  hhif  dbtilfliU  whih  ill 

state  variaoles  except  the  temperature  are  held  fixed.  The 
first  part  of  thts  e»pi>*Ns:on  follows  from  me  internal  energy 
Iprm  -en:!-  the  cecen.d  tors  cfsial  to  the  ent.ropy  ot  the  system 
the  if.  teyrai  of  the  neat  iiiput  over  the  temperature.  The 
.htegratior  :s  performed  over  the  volume  of  System  I.  Equa- 
tio*”  (14)  4lso  be  written  m  the  form 


V.  ..  j  [f® 
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If  9  <<  7^,  that  iS,  small  departures  from  the  equi  I  ibrlpm 
state,  V_  cap  be  wiittep.  as 

'/c  =  i  i  av  (It,) 

r 

Therefore,  the  function  V  can  be  written  in  the  for® 

^  *  2  i  dv 

r 

where  is  the  value  of  V  for  ar  isothermal  process  for  example, 

IP  thermo  elasticity  the  value  of  V  will  be  the  value  of  the 

r 

isothert.al  free  energy  or  strain  energy  integbatcd  over  the  volume 
and  c  Will  be  the  value  of  the  specific  heat  at  constant  stislnu 
In  this  wav,  the  effects  of  the  temperature  and  the  remaining 
state  variables  Ct?'  ce  sidered  separately 

Now  that  the  s i gn  i  f i car  •  e  of  the  thermal  potential  has 
beer  dismissed,  the  introductior  of  external  disturbing  tbfees 
can  be  corsidereg  m  order  to  formulate  the  complet*^  Lagrangian 
aguar.oi  1  Ml  T.!  e  slate  variables  in  the  most  gt^neral  case 

rhermoOvram,  Sy'.^em;  t-ifernal  lories  Applied 

’■pe  eguatior  for  the  time  histories  of  *he  state  . 
vaiiaCies  of  t^e  t  hermod  yr  am  i  *  ivstem  were  giver  for  th#  ,'c'asu 
or  sc.-c  nyf-'r-ai  forces  by  Egs  (H).  In  order  to  bbtim  tha  u" 

f  rm  of  tre  eguaticrs  for  the  case  where  extomal  forcet  a,r.b 
a-ti.g  on  Sv-Ttetn  i,  the  entropv  of  the  total  sy  stem  w, ’»  I  again  ' 
be  expres  ed  in  terms  of  the  entropy  and  ifie  internal  entf.gy 
cf  System  '  ’’ms  will  parallel  •  he  d  l  scu  s  *i  or  of  the  •mq- 

'  M  ^  ■  ' 

nitiia*';#  of  fhe  ^r-ermal  potential  fupvtlor  /  However,  m 


"If**?' 


this  case  conservation  of  eneigv  will  introduce  the  effects  of 
the  external  torces,  Cons-.der  '  heat  reservoir  at  temperature 
-  9  adjoining  the  primar/  System  I  and  let  the  heat  flow  from 


this  reservoi’-  be  equal  to  Conserv 

3  *s  then  expressed  bv 


ation  of  energy  for  System 


u  =  h„  ■  h 


<*'iere  h  is  aga-r  neat  flow  from  System  1  to  the  System  II,. 

The  er-troDV  S'  of  the  total  svstem  is 


S  =  S  -  S,,  -  S 


II  'III 


n'-ere  S,,,,  rs  the  ertropv  oi  the  reservoir  at  the  temperature 
^  '  "'"■i  equaf.o-'  ca"  be  written  m  the  form 


S  -  ?  •  fh,/T^  r  0.)  ♦  (h/T  ) 

*  *  r 


w'  :ch  t  eiomes  '.  pon  » !■  1 1 ocl  . c 1 1 or  of  Eg  *  16) 


S  ,  _  !  '  •  h  'f  I  /- J  li/r  „  Q  j  ] 


If  i'’  <<  ’■  ,■■■  it  •» 

^  V  ' 

s  V. .  •  e* 


e*-  -Oilcw,  that  fhw  entropy  of  the  total 
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3'  =  S  -  (U/T^)  ^(hje/T^) 


The  last  tei»  ip  this  enpressio''  *ili  be  behoted  as 
where  the  ter® 


Will  be  called  the  e''tTcpv  inflow  to  System  I  Therefore* 
the  expression  for  tne  entropy  S’  of  the  total  system  can  be 


wr 1 tter  ir  the  ♦orm 


’3=  •  V  ♦  S.O 
r  ? 


VIT) 


where 

V  =  U  -  T  S 
r 

’r-e  rerm  S,,0  is  .la.-lai  fo  a  generali^ecll  work  e*- 
pressiO''  in  the  '.aarangsan  tormulatioh  of  the  governing  eqt-a- 
‘•nrs  i  lie  ch  j<r  . 4 1  >vste*  The  quentity  Sj  is  the  force 

•--ct.on  and  0  •  ■  'he  vonti^gate  displacement  3  milar  e»pre»- 
sions  for  rh*  xiors  quantities  nf  a  thermal  system  in  terms  of 
phys.^Hi  forces,  cohcen t r a t i pn  giadients,  wc;  ,  can  be  deter’ 
m.ned  aoj  some  t*f  these  add  >  1 1  o' a  I  cases  are  outlined  in  Pet  t  '  i 
In  view  yf  tn ;  ,  .r  i*  poss.b.e  to  write  hq  M ' '  m  tie  general 


•  orm 
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r,s'  =  -  V  .  0,,, 

.wnere  is  the  generalized  force  conjugate  to  the  generalized 
displaceaent  For  exanple^  the  generalized  force  conjugate 

to  the  tenperaturet  as  outllneu  above^  is  the  entropy  flowc 
Onsager's  relations  can  be  applied  to  the  complete 
system  consisting  of  the  primary  system,  the  heat  reservoir  at 
the  equilibrium  tes'perature  and  the.  energy  reservoirs  equi¬ 
valent  to  the  applied  external  force:.  This  application  is 
similar  to  the  case  where  no  applied  femes  are  acting  on  the 
system,,  In  view  of  Eqs  (4l  and  (6),  the  equations  for  the 
state  variables  take  tne  form 

T  51_  ,  4fi. 

r  ^  ^  ^  q 

or 


Ji± 


Equations  Msl  are  ♦h*  governing  equations  for  the  state 


(18) 


variables  q^  of  the  therma>  svstem  These  equations  are,  to 
repeat  the  statements  made  above,  similar  to  the  equations  of 


a  mechanical  system  in  which  V  is  the  potential  energy,  0  it 
the  dissipetiun  function  and  the  ere  now  the  extornal  for¬ 
ces  act  ing  on  the  svstem  Ac  additional  way  of  writing  Eqs'  (18)  is 


-  -jL  f-i  dill 
^q^  *■  ^  '  2  dt  ^ 


wh'.c*-  expresses  forces,  both  txtetrel  ano  internal,  m  terms 
of  the  entropy  production  m  the  system  Adoitiofal  properties 
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of  Eqs.  (18)  are  also  discussed  in  Refs-  [2,  3,  ?]  but  these 
will  not  be  discussed  here-  Insteaui  the  equivalence  of  the 
equations  to  a  variational  principle  will  be  shown  with  a  view 
toward  obtaining  a  variational  principle  for  heat  conduction 
analysis. 

The  function  D  can  be  written  in  the  following 
operational  form 

D*  =  ( 1 /2)  pb^^  q^q^ 

where 

p  -  d  ''d  t 

If  the  d.ssipjtion  function  is  written  in  this  form,  cqs.  (18) 
ca''  b“  written  ir  the  variational  form 

♦  50*  =  Qj  5q^  (19) 

variation  i»  with  respect  to  the  gtntrelixed  coordinetts 
ar,(  t>#  , operator  p  :*  treated  as  a  constant  when  calcu‘> 
iai.  ’g  me  var  ation  'his  form  of  the  varistional  principle 
IS  convenient  m  some  cases*  but  the  varietional  principle  can 
t# , «s tab  I i sned  also  on  the  easi«  o*  the  definitions  of  the 
dissipation  function  eng  thermal  potential*  Ref  td)  T^is 
Will  be  the  method  ter  tn#  determin«tlon  of  the  form  of  tho 
va’iats»n«i  principle  ir.  heat  coRductior  analysis  to  b* 


disttisseo  "e«t 


Flow 
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The  Lagrangian  equations  obtained  above  are  equiva¬ 
lent  to  the  variational  principle  o;ven  by  Eq.  (19).  Instead 
of  continuing  the  study  of  a  theraodynaam.  ^ystea  in  coaplete 
generality,  the  particular  case  for  heat  flow  alone  will  be 
studied  ti^it  IS,  the  eiain  concern  Wi'l  be  with  the  theraal 

field  0  In  particular,  the  expressions  for  the  dissipation 

function  and  thermal  potential  will  be  found 

’’ha  r-eraal  puiertial  for  this  restricted  case  fol¬ 
lows  lesiediateiy  troa  cq.  «  Ih' .  i^d  is 

V  s  t  dv 
•  ■  '  r 

’*'e  3:  .oatiwin  t.. -ct  ior  *■  ;>  e»prtsied  m  the  form 

I  -  \  ‘£'^  '  X  ?Ha*e  nf  f-tropy  f* reu  'itton) 

•herefore  it  .x  netes^arv  to  express  the  rate  or  entropy  pto- 
■iv  ir  teres  n»  Vnci,r  q!»antiti*e»  of  the  thermel  tyetem. 

Th , »  van  b«  dc'-e  t»  applvine  first  and  second  laws  of  ther- 
etuiyn  ^  to  v.-h  4  syvte*  ’“he  first  lam  i* 

i2 

dt 


where  U  ib  the  internal  energy  of  the  system  and  H  is  the  heat 
flux  vector  field,  that  is,  the  heat  flow  per  unit  tine.  The 
second  law  gives 


dS  _  dU 
d  t  ~  d  t 


div  H 


(20) 


wh*re  T  i£  the  lenperature  and  S  is  the  entropy.  Equation  (20) 
car  be  written  m  the  form 


^  -  div  (H/T)  =  -  (H/T^)  VT  (21) 

’’be  left  hard  side  of  this  equation  is  interpreted  as  the  change 
;r  entropy  of  a  unit  volune  plus  the  entropy  flow  across  the  sur- 
*sc.e  of  this  volume  Thv  ret  change  in  entropy  for  the  unit  volume 
IS  then  given  by  the  right  hand  side  ot  Eq*  (21)*  this  term  is  the 
entropy  produrtmo  in  the  volume.  furthermore  if  the  thermal  force 
VT,  s  linearly  related  to  the  flux  field  such  that 

H  *  -  VT 

tne  entropy  prodvution  term  proportional  to  the  square  of 
the  heat  flu*  field  Tne  minussigr  the  above  evpresr.lon 


arises  since  the  heat  flow  must  be  opposite  to  the  temperature 


gradient  Therefore,  the  dissipation  function  can  be  expressed 

in  terms  of  the  squire  of  the  neat  tlux  field 

The  remaining  term  discussed  with  reference  to  the 

general  thermodynamic  system  was  the  generalized  force 

This  force,  in  the  case  c*  •?  purely  thermal  field)  is  immediately 

‘"'own  from  tne  previous  discijssions  and  it  is  equal  to  the  heat 

flow  field  Oiv  cse;!  by  ’’  ihe  factor  T  used  above  can  be  drop- 

r  r 

ped  since  it  is  common  to  all  terms  in  Eqs  '18) 

■’'he  ^^nctiors  m  tne  Lagrangian  equations  can  now 
be  e^precsed  m  terms  of  the  temperature  field  and  the  time 
rate  of  change  of  the  fieat  flow  field  Once  so  expressed)  as 
was  do''e  in  the  above  paragraphs,  the  discussion  for  the  general 
tn • rmodvr am ; ;  system  indicates  tne  form  of  a  variational  principle. 
wh'.cn  :>  applicable  to  a  thermal  system,  m  particular)  to  the 
t-  H\  0*  :  on  -  •*  "' *■  ....'•«••>  Tn.s  variational  piin- 

*  pie  was  o’.^'e'"  hv  Rofo  f5)  6,  and  it  was  shown  to 

be  equival«-(  ♦  to  tne  heat  conduction  equation  Once  this  equl 
vilence  was  established,  the  heat  flow  field  related  to  the  tem- 
peratoiP  f,eld  hy  conservation  ot  en'-rgy  was  expressed  in  terms 
of  generalized  coordinates  ''’he  variational  principle  was  then 
shown  to  lead  to  Lagrangian  equations  for  the  generalized  co¬ 
ordinates  of  tne  temperature  field  In  this  *»ay  the  motivation 
behind  the  beat  conduction  variational  principle  as  being  related 


to  th*»  Lagrangian  representation  of  a  general  thermodynanic 
system  is  shown 

Further  discussions  of  the  principle  for  heat  con¬ 
duction  were  giver,  m  Ref,  [13]  together  with  a  number  of  ap 
plications.  A  further  application  ui  the  principle  will  be 
giver,  isi  the  next  section.. 


SECTION  V 


THE  HEATING  OF  SLABS  EXPOSED  TO  TIME -DEP ENDENT 
HEAT  I  LUXES 


S'? 


The  problem  of  the  determination  o^  the  temperature 
distribution  in  one-dimer. sional  slabs  heated  bv  arbitrary  time- 
dependent  heat  fluxes  has  received  much  attention  in  the  last 
few  years.  This  interest  has  been  er gendered  by  problems  aris¬ 
ing  in  the  design  of  structural  components  tor  re-entry  vehicln^ 
In  these  design  situatiors.  the  .iface  temperature  is  small 
compared. to  the  stagnation  temperatuie  ot  the  air  and  thus  cuo 
be  neglected  m  the  usual  aerodynamic  ►eating  rate  This  as¬ 
sumption  then  reduces  the  neating  rate  at  the  surface  to  an 
,arh  1 1 r a-r y  functioh  of  time 

The  e  *ac't .  ■so  I  j  t  ibn  to  th*  problem  of  th#  heating 
0  ‘  f  in  i  te  s  I  abs  and  semi  •  in.'i  •,  n  \  te  Bodies  bv  an  arbitrary 
t  ime  "dependif^i  t  heat  t  lu.k  j  s  .known  j.. Ret  $  [40  to  bO],  and  this 

tolutioc  '’an  5e  expres'i^d.  ;r- '  terms  of  ;  . the  solution  tor  a  ton- 
sti)"’  '■♦at  f  i Ih'e  ^special  case  o»  apaiaPoi;.  neat  inp>  t 
wa-i  treated  by  Ks'iie  ,;}ha  Ye.n,  pr*  ('ulll  while  the' general  '.!»%• 
or  tol'.  i  sal  protiles  w., ..  .  t  reate  ihih.- Rpf  [h*!^,:'h?j  Si -jO  ;  an 

and  t't  1  angu.i  at  hee't  pul  v’e«  were  tnyesfigatedi’n  Ref.;  |b4} 

..THe-s*  io  i  a  t  •  were  .  t  a  1  nft  .  by  sr.vd''r(j  tn*  >>'euat  rressscal 

h  ea  t  ■von<j'..f' t  i  un  i^r.nb  i  #«■-;  and  a*#  obtained  r'-  a  ’  ;  r'tv.  tn'at  ;  s 

V.':.  li.'-  :  n;.-u- <  n  t'  f  0  r  hufeerital  eva  I  ;•  a  t  >  si” '  A*.  ,  a- 'con  segu  en  » 


3b 


ar  approximate  method  solution  is  useful  for  obtaining  the 
required  temperatu-o  strib*  tio'' s  in  an  expedient  way. 

The  method  to  be  employed  here  is  8iot*s  variational 
principle  for  heat  conduction  Tnis  principle  and  its  appll' 
cation  to  problems  with  neat  flux  boundary  conditions  was  dls- 
.uased  in  detail  ir  Ret  tl3]  The  present  discussion  of  the 
principle  applied  to  tne  problem  of  the  heating  of  slabs  with 
t.me  dependert  fi(.»es  will  be  based  upon  Ref.  Cl3j-  It 

*11  b*  ass.meu  ttat  the  reader  is  familiar  with  the  method 
or  soiut:o;  Outlined  m  Ref.  [13]. 

problem  will  be  formulated  for  an  arbitrary  heat 
.'••pcx  and  then  p  ar  ti  cu  I  a  r  :  z  ed  to  the  case  of  a  triangular  heat 
p'il'ie  As  was  discwssed  ir  Ref  Cidlx  the  formulation  of  the 
proD*em  for  the  heating  or  slabs  proceeds  in  two  phases. 

'  .  I  f  1  r  :  te  Solid 

■'•■e  first  phase  of  heating  coiiespohuS  to  the  semi- 
portion  ''‘r.  i  s  phase  ends  when  the  penetration  depth 
s  equal  to  the  slab  tnictfness  The  temperature  distribution 
;s  assumed  to  parabolt:,  'ig  I  .and  expressed  in  the  form 


0  q  f  1  -2-1  ’ 


I  • 

■•1 


■  e  I  e 


xs  the  Surface  temperature  and  Qj  is  the  penetration 


I 


depth-  The  generalized  coordinate  q2  is  selected  as  the  in¬ 
dependent  coordinate  related  to  through  the  expression  for 
overall  energy  balance,  Ref  [13].  Evaluation  of  the  appropriate 
derivatives  for  the  thermal  potential  and  dissipation  function, 
together  with  the  thermal  force  and  surface  heat  flux  leads  to 
the  following  exor-issions  ■. 


=  (  C^/k)q^  q^Ci  l/42)q^q2  +  {  1 3/315 )  q  ^  q2  ] 
Q2  -  cqJ/lC 


--  (  c/3)(  q^q2  r  q,  q^) 


Introduction  of  the  flu«  cot"*ltion 


=  hql  t) 


wneic  q(ti  is  a  d rmen •  j on  1 es t  function  of  ti«»,  and  the  para¬ 
meters. 


n  ■  5  2’--  ® 


‘■t  =  ijk/ho 


T  •  at/b 
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will  reduce  the  governing  variational  equation  and  flux  condi¬ 
tion  to  the  following  dlaensionless  equations! 

*  2 

15i^f|  +  =  147\|»  (22) 

d(  ♦l^).^dT  =  3q(T)  (  23) 

Equation  (23)  can  be  integrated  iamedlately  for  the  Initial 
conditions  f(  =  >^  =  0,  atT  =  0  and  the  function  elialnated 
between  Eqs-  (22)  and  (23).  The  resulting  equation  is 

Hnn  i^>n^Co(-)/9(t)  -  g(o)3  s  147  (24) 

where 

gi^r)  s  f  q(T)dT  (25) 

end  9(0/  is  the  value  of  g  at  t  *  0.  If  the  eubstitutlnns 


i 


*1 


2 


r  (t)  *  q(T)/9(T)  -  9(0) 


(26) 


41. 


are  Bade,  Eq.  (24)  will  reduce  to  the  equation 

2  +  (30/il)2F(T:)  =  (  294/11)  (27) 

Ihe  solution  to  Eq.  (27)  under  the  initial  condition  z  =  0  at 
T  =  0  is 

2 

2  =  q  =  b  Jexp[ajF(T:)dT]dT/exp[ajF(T)dT] .  (28) 

0 

where  a  =  (30/ll)  and  b  =  (294/11). 

This  solution  gives  the  penetration  depth  as  a  function 
of  time  for  arbitrary  heating  rates.  Once  q  =  q(T)  is  found, 
the  surface  temperature  history  can  be  found  from  Eq.  (23). 

The  resull  can  be  written  in  the  form 

i>(t)  s  3Cg(T)  -  g(0)]/q  (29) 

If  the  heating  rate  is  such  that  g  (0)  equals  *ero, 
for  <?xa«pie,  a  polynomial  function,  then  Cv,.  (28)  reduces  to 
the  torm 

C9(t)]*  dT/Cg(-> )  ]• 

0 


(30) 


A2. 


In  this  case  the  solutions  for  the  penetxaLion  depth  and  surface 
temperature,  Eqs.  (28)  and  (29),  assume  a  simpler  form.  A  parti¬ 
cular  form  of  interest  for  the  heat  input  is  the  case  where  the 
heat  input  is  of  the  form 

q(T)  ~  Ct” 

that  is,  a  heat  flux  dependent  upon  a  power  of  the  time.  Evalua- 
2 

tion  of  ri  from  Eq.  (30)  yields 

f|^  =  294-1/41  e  30n  (31) 

This'result  shows  that  the  penetration  depth  is  independent  of 
the  constant  C  In  the  heat  input  and  varies  as  the  square  root 
of  time.  The  factor  of  proportionality  Is  dependent  on  the  power 
of  the  trim*  of  the  heat  flux.  Equation  (31)  is  an  interesting 
result  that  compares  with  the  result  obtained  by  diet  in  Ret.  [9]. 

In  this  reference,  Riot  found  that  if  the  surface 
lemperatvre  varied  as  -r^',  the  penetration  depth  relation  i# 

*1^  =  147t/13  Ibn  (32) 

These  two  n^sui  ts,Eqs.  (31)  and  (.32),  show  that  the  p 'jnetration 
depth  varies  as  the  -inu-jre  root  of  time  for  two  casest  '.urfac# 
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temperature  dependent  on  a  power  of  time  and  heat  flux  dependent 
on  a  power  of  time. 

If  n  =  0  in  Eq.  (31),  that  is,  a  constant  heat  flux 
and  n  =  (l/2)  in  Eq.  (32), also  a  constant  flux,  the  penetration 
depth  relation  is 


=  294t/41 


This  result  r'or  constant  heat  flux  was  also  noted  in  R»f.  [13]. 
If  n  =  1,  a  linear  heat  input 


q{i )  =  t/ 


where  equals  a  constant,  the  penetration  depth  relation  is 


tl‘  =  2v4t/71 


The  rnrrespondvnq  surface  terptratuVe  history  from  Eq..  (29)  is 


)  >  (  3/2Tj)(  7  1^294) 


L^ustlu.i  {-33)  corresponds  t&  the  flr^t  ot  *  triangular  heat 
pulse  and  Eqs,  { 34)  and  ( 33)  are  the  corresponding  .penetrallvn 
depth  and  surface  temperature  histones  iheh  e  s  -i,  he  ■ 
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penetration  depth  equals  the  slab  thiuKiiess  b.  The  penetration 
time  from  Eq.  (34)  is 

Tp  =  71/294  (36) 

and  the  surface  temperature  at  the  penetration  time  is 

\|)(Tp)  ^  (3/2Tj_){Tp)^  (37) 

Equation  (3/;  gives  the  surface  temperature  at  the  penetration 
time  which  is  to  be  used  for  the  initial  condition  for  phase  2, 
the  finite  sleb. 

Finite  Slab 

The  temperature  distribution  for  the  second  phase 
corresponding  to  the  finite  slab  is  shown  in  Fig.  1  and  is 
written  in  the  forn 

©  =  {  Q  .ni  , I  1  -  ^ ^  *  q  . 

i  ^  b  ‘4 

If  the  analysis  given  in  Ref.  [13]  for  finite  slabs  is  applied 
together  with  the  paremetert 


/ 


-15. 

T  =  at/b*" 

’1.3  =  ^1,3 

the  variational  equation  and  flux  condition  become 

•  * 

lOm^  +  32^2  =  -  V3) 

•  • 

-*  2^2  =  3q(T)  (38) 

The  flux  equation,  the  second  of  Eqs.  (3B),  can  be 
integrated,  and  the  function  eliminated  between  the  equa¬ 
tions  to  find 

e 

’1  *  ^^’1  “  ^  ?iCg(T)  -  9('fp)]  * 

Equation  (39)  can  be  inteqrated  to  yielo 

♦  j'‘»i(^p^f«*pUlT)(8q(T>2lCg(T;-g(tp)]f?*j(Tp)}dT/exp{21i)  (40) 

This  result  gives  the. surface  teapdrature  as  a  function  of  time 
for  an  arbitrary  heat  jflput.  The  rear  temperature  history  can 
be  found  immediately  from  the  flua  conditloh.  If  the  flux  is  a 
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fiolynoalal}  tor  axample)  Eq.  (40)  can  be  integrated  in  closed 
fora  for  the  surface  temperature. 

In  suaaary  ther^  tb.e  basic  equations  for  an  arbitrary 
heat  input  are  given  by  Eqs.  (23)  and  (29)  for  the  semi-infinite 
solid  and  by  Eq.  (40)  for  the  finite  solid.  These  equations 
shon  that  the  solutions  for  an  arbitrary  heat  input  can  be  ob¬ 
tained  in  closed  form  for  many  cases.  An  example  of  the 
use  of  this  method  mill  be  indicated  to  show  how  the  surface 
temper a tu r e  h i s t nr i es  are  found 

The  example  to  be  treated  is  the  heating  of  a  finite 
slab  by  a  triangular  heat  pulse.  This  example  was  treated  in 
Ref.  Cbdj., 

Triangular  Heat  i'ulse 

The  heat  pulse  to  be  considered  it 

»  2  “  *  T  ia  (41) 

\  -  * 

The  s^urface  temperatures  of  a  finite  slab  of  thickness  b  are 
required  for  r^a.  ihe  results  for  the  s*mi<^inf  Ini  te  portion 
are  given  by  Eqs.  (34)  to  (3T)  for  *  1.  The  penetration 
timt.Jftj?*  Eq.  (36)  is  • 
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1  ^  7l/2‘'-A 

P 


This  time  corresponds  tc  the  starting  time  for  the  second  phase. 
Equation  (39)  in  this  case  reduces  to 


21*^  = 


10.  5t 


T>T. 


Th  e  solution  t  s 


T*  •  T  *  ^  -  I  CCS6e»p{-21't) 


( 42) 


rom  #hich  tne  inner  surface  temperature  is  found' 


2t3  . 


*3  *  s’*'  *  6  ^  *1^  i028e»p(  «21t  )  (43) 

E  juit  s  jnv  42;  («3'  ate  the  surface  temperature 

histones  These  solution^  rpply  until  -  *  at  mhlch  time 

the  fiu«  condition  iaaes  •  diiferent  form  from  Eg<  (41).  At 


this  time 


qi't)  -  2  - 


and  tne  gove-m-nq  eq.iatiur.s,  Eqs  (38),  are  again  integrated 
for  tne  initial  conditions  corresponding  to  t  =  1  These  c'n- 
ditiors  are  evaluated  from  Eqs  (42)  and  (43)  at  t  =  1.  Evelua 
tion  leads  to  the  values 


0  S  1  46 


1  ;  =  0  341  2'7 


to  t'j  used  as  the  ti  al  cord,  t  ions  Solution  oi  Eqs,  (38) 
'■■(‘'r  these  ir;  t'.al  conditions,  vields 


.  1  r, 


'3rU6etpr.2l(-  :)]  (44^ 


■■1  f 


0  Oldb ) j  (  4h) 


tqn.)r^-''is  !  44  '  and  ,  4h )  q-.e  me  urlace  temptri  tures  of  tfu 


^  '  a  0  *  0  .  ?  ‘  e  t  I  iTi  e  ’  ,i "  q  - 


eh;  I#  Fqt;  (4y'.  joJ  (43)  give 


the  s..!ia.,e  ’  ••tip  e  t  a  t  m  e  .  m  r  th-  ti*.*  t  -  s  ..  th* 

"  •  -  -  p 


inner  surface  temperatute  i.c  of  co  rse  equal  to  zero  until  t-e 
penetration  time 

The  above  ies.;it;  for  the  surface  temperature  n.  c  shov. 
in  Fig  2  tog  c  tn  1  with  the  r  es;.  1  fi  f  rom  Sef  [54]  It  is  -  een 
from  th*5  tiqure  t.-.ir  the  present  methoc  is  accurate  as  an  ap¬ 
proximate  (netf'o.J  ot  solution  rot  the  case  o?  a  triangular  heat 
pulse  ■  . 


50 


SECTION  VI  CONCLUDING  REMARKS 

A  discussion  and  application  of  Biot's  variational 
principle  for  heat  conduction  Mas  presented  in  a  previous 
paper*  Ref,  Cl3]  This  reference  contained  a  dlsccssior.  of 
the  principle  from  a  mathematical  viewpoint  together  with  a 
discussion,  of  tne  method  of  solution  using  the  principle 
A  number  of  different  one-dxmenaional  heat  condurtlon  problems 
were  treated  to  show  the  applicability  of  the  method. 

In  the  present  paper*  the  literature  discussing 
Biot's  variational  principle  is  reviewed*  in  addition  to 
other  variational  principles  and  approximate  methods  of  solu> 
tlon  for  heat  conduction  problems  The  thermodynamic  founda- 
tlons  of  the  variational  principle  are  reviewed  in  Section  IV 
showing  that  the  variational  principle  for  heat  conduction  has 
Its  basis  lo  the  physical  concepts  used  in  irreversible  ther¬ 
modynamics  The  latter  part  of  Section  IV  then  shows  how  the 

form  of  the  principle  for  heat  conduction  follows  from  the 
general  formulation 

An  application  of  the  method  employing  the  Ideas 
given  in  Ref.  [13]  is  presented  in  Section  V,  This  application 
involves  the  heating  of  slabs  exposed  to  time-dependent  heat 
fluxes  For  this  type  of  heat  flux,  the  present  method  of  solu 
tlon  admits,  in  many  cases*  a  closed  form  solution  for  the  sur¬ 
face  temperatures  The  case  of  a  triangular  heat  pulse  Is  then 
presented  In  detail  to  show  how  the  method  is  applied  for  pro¬ 
blems  of  this  type 
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FIG.  I  TEMPERATURE  PROFILES  IN  TERMS  OF 
GENERALIZED  COORDINATES ;  HEAT  INPUT 
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FIG  2  SURFACE  TEMPERATURE  HISTORIES 
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Wilght-Patterson  af  lta<» 
Oh  1  n 


ARl  (afro) 

ATTNt  RR.,A  —  Library 

Wrlght-Patterson  AF  Base 

Oh  1  o  (2) 

AFIT  (AU; 

ATTNt  MCL  1  ••  Library 
Wright-F.tterson  AF  Bate 

Ohio  ^  ^  j 

APGC 

ATTN.  PGTRl  -  Tech.  Library 
Eg)  in  AF  Base 

Flo'^lda  (Ij 

AFCRI.  (sFRD) 

ATTNt  :rre  (2) 

L.G.  Mantrnmb  Field 
Uedo'rt^  Nesssuhuielts 

P.  Box  AA. 

Wrlght-Pattersen  AF  Base 

(1) 

AMI 

Dt  rec  tor 
BM 

AT  IN  t  Library 

Aberdeen  Proving  Ground 

Nar y land  ( j  j 

Office  'tf  Ordnance  Research 
Box  CM 

Dube  Station 

Durham,  North  Carolina  (1) 

Army  Rocket  &  Guided  Missile  Ag«n 
ATTNt  Technical  I Ibrary 
Redstone  Arsenal 
A 1  abana  ^  ; ) 

'tlgnal  Corps  tnqlneertng  lab. 
•mTNi  3  IGi^M  1  f  l-RPU 
'  r  t  Monmou  th 

Newfarsey  ,jj 


OffiCfi  of  tne  Chief  of  R&D 
ATTM;  Scientific  Information 
Deoaftrrent  of  the  Army 
Washincton  25,  D.  C.  (] 


Office  of  Naval  Research 
ATTN  I  Mechanics  Branch  (1) 
Airbranch  f’l) 
Washington  25 ,  D.  C. 

Naval  Research  Lac.-c  t;  ry 
AlINj  Documents  Library 
Washington  25,  D.  C- 

Naval  Oic'nance  Lafcoratuiy 
A  T  TN :  Library 
White  Oak,  Maryland 

David  r -i  y  '  !  Mode.  Pa c  ;  n 

Aerodynamics  Laboratory 
A  T  T  N  1  i  b  '  a  r  y 

Washington  "7,  C.  C. 

Chief,  Bureau  of  Ordnance 
Depa *  imen t  of  the  Nacy 
ATT'>:  Special  t  rejects 

Office.  SP-2:22 
W  a  s  h  i  n  g  t  j  M  25,  D .  C . 

CldiB  nCVCRNMLNT 


A  T  T  N  1  D  o  c  u  in  e  a  t  I  I  o  r  a  t  y 

15  20  H  ')  ♦  r  "  \  .  A 

Wasnlitgton  25,  C.  (5) 

Nations;  Swr'.cu  '>tandaius 
U.  5.  Depar’ment  of  C  ammo 
AiiNi  Lin^ui  ■ 

Washington  25,  T.  C.  (1) 

Office  of  Tosh nica!  Servlf's 
'J.  S.  Depaifment  of  Commerce 
ATTSi  Te..nnical  Reports  'ierllon 
Washington  25,  0.  C.  (  11 

Natlo.nal  science  rounOation 
ATTNi  Eng.'g.  Scl.  Dlv, 

1951  Constitution  Ave.  ,N.W. 
Waihl.ngtor.  25,  0.  C.  :  I ' 


(J.  S.  Atomic  Energy  Corri mission 
Technical  information  S-- 
1901  C  0  n  1 1  t  u  t  i  0  n  A  ,  N  .  w  , 
Washington  25,  D.  C.  (l) 

U.  S.  Atomic  Energy  uomrri.  ■  'n 
Tei.I.nic.il  In  i  u: -sy  il  an  L  x  v.  eu  s  i  on 
P  .  0  Box  6  ?. 

Oak  Ridg«i,  Tenness.  e  (1) 

JOURNALS 

South'eost  Research  Institute 
ATTNi  Applied  Mechanics  c  .‘k 
8500  Culebra  Road 
San  Anton'o  6,  Texas  >  .' ) 

A  e  1  o  n  a  u  I  i  f.  a  !  F  n  a  i  n  e  n »  r,  g  Review 
2  East  b^th  Street 
N  ow  Y  0  r  k  2 1  ,  N  ,  Y  .  (  1  ! 

Inst,  of  A  e  r  o  a  a  u  t  i  c  a  1  c  c  i  i  .  s 


ATTNs 


i  ^  r  a  r  y 


2  fa’'*  b^th  Street 
N  e  w  Y  n  r  k  2  1  ,  N  .  Y  .  (  1  ) 

^  Oht:  IGd  jY R I ;  AN  t  T  ;;  ny  ■ 

Chairman 

Canvilan  Joint  Sta'f  '  DRP/r";  1  S  ' 
24*. 0  Massachuset  ts  Ave.,  N.  A. 

W  a •  n  I  n  g  ton  25  ,  D.  C •  t  1 ) 

5  I  ’  p  r  L 1 1 1 

I’iatlunai  Aeronautic. il  r 

C  t  t  .1  e  a  ,  O')  !  r  I  *  n 

C  a  n  a  d  a  1  ) 

Unlv.rslty  of  L'rnnto 
institute  n  f  y  "  r  0  p  h  >  •  I  c  s 
Alt’)  I  library 

I o  r  o  n  t  o  5 ,  Canada  (  I  ) 

Cen.  for  Ixper.  Aeiody 
AlTNi  Libtery 

Hho  I  e  -  Sa  i  n  t  -uene  b  »•  (  Hr  1  q  1  ._,u  »  i 
72,  Chaussee  de  W.ler..,. 
Brussels,  Belgium  (l' 
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educational  INST  1  T'JTIOMS 

Auburn  Unlvpr*^ltv 
Departnenl  of  Mer.h.  Lnqiq. 
Auburt),  Alabama  (1) 

Brown  Unlvrrslty 

Gifts  and  Exchanqi?s  lltrary 

ProvJoem.*?  12.  Nho-Je  I- land  (1) 

University  of  California 
Englneerlnq  Oepa.tment 
ATTN  I  Library 

Los  Anneles  24,  Califotnia  (ii 

California  Inst,  of  lerh. 

ATTNj  Jl’L  Library 

4800  Oak  Grove  Drive 

rdsadena  4,  Cc’.lfurnla  (1) 

California  Ini...  of  Tech. 
Guqgenhelm  Aero.  lab. 

ATINi  Aeio  lie.  (Haute  to 
Prof.  M.A.  Liepmann) 

!’as adona  4,  California  (1) 

Colo:  ado  ..tatt  ’Jr.  1  v  e  r  s  i  t , 

“  p  t  ■  of  t.  <  V  1  1  t  n  g Q , 

ATFNi  Prof.  J.E.  Cernak, 

a:ce 

For  t  Co  1 1  ln»  ,  Col  oral.*  (l) 

Columbia  University 
Dept,  of  Cl<,  Etiq'q,  A 
Engrg.  SSech* 

AITNi  Lib.  (Route  to 

F  t  o  f  .  "  .  "  e  ' 

New  York  27,  K.  y.  (  i ) 

University  of  Florida 
rniine"rlr.g  Mem.  Dept. 

ATTNi  I.  i « !  d  r  Y 

'jaln«vl;le,  Florida  (i; 

Ha rva id  Un 1 vei > i i) 
eept.  of  Tnorq.  Si  If 
ATlNi  Library 

Caabrldqe  38,  Mas- .  (  • l 


John  Crer,>r  I  ihrary 

8&  E.  Randolph  Street 

C(',  icaqo  1,  Illinois  (j: 

The  Jo.brjs  Hopxlns  Ur.lveijlty 
tppHed  Pnysir  5,  Laboratory  lib. 

96  21  Ge(jr.,ij  Avenue 

Si Iver  Spring,  Maryland  f  1) 

The  Johns  HepPms  University 

Department  of  Merhanlrs 

ATTN:  Lioraiy  (Route  to  Profs, 

!  .  ' .  C 1  o  0 ..  e  1  &  3 .  C  >■  r  r  s  i  n  ) 
Ba.'  ti.sore  18,  Maryland  (  i) 

Massachusetts  Inst,  of  Tech. 

ATTN:  library  (Route  to  Dept, 

of  Aero  &  Mech.  Engrq. 
and  Dept,  of  Anp.Mech.) 
Cambridge  39 ,  Mass.  {  I  j 

MMw'st  Restarch  Institute 

ATfNi  Library 

42I>  Volker  boulevard 

►  mr  as  City  10,  u:.,sourl  (  '  ) 

North  Carolina  State  Colltje 
Div,.i,,n  jf  tngrg.  Research 
ATlRt  Technical  Library 
Ral-igh,  North  Carolina  (  1) 

t,..ytechnlt  Institute  •  Bklyn. 

<  I  l  N  :  i  l  b  I  a  I  y 

113  Jay  Street 

o  1  y  n  ; ,  Ii ,  Y , 

'•nnsylvanla  Stale  Univeralty 
Dept,  of  Aero.  Enqlnetrlnq 
*iIN«  Ltbraiy 

tn I ve  r  s I ty  Park  ,  Pa.  ( 1 ) 

ihe  lames  Eoirestal  Res.  Canter 
>  •  1  n  •  t  on  fn  ;  .e  »  si  ty 
'  I  " . ;  !  •  b .  (Route  to  Prof. 

Bogdonoff) 

i  :  1 1.  i  *  >o  ,  N  .  j  .  i  :  ) 
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Princeton  Univtrslty 
Dept,  of  Aero.  Engrg. 

ATTMi  Library 

Princeton ,  N.  J .  ( 1 ) 

Rensselaer  Polytechnic  Inst. 
OepartJient  of  Aero.  Engrg. 

ATTNt  Library 

Troy,  N.  Y.  (  1) 

Stanford  Research  Institute 

Documents  Center 

(at IN I  Acquisitions) 

Menlo  Park}  California  (1) 

Stanford  University 

Dept,  of  Aeronautical  Engrg. 

ATINi  Library 

Stanford}  California  (1) 

Defense  Research  Laboratory 
University  of  Texas 
P.  0.  aox  8029 

Austin  12,  Texas  (1) 

New  York  University 
Institute  of  Math.  Scl. 

ATTN:  Library 

New  York  3,  N.  Y.  (1) 

Yale  University 
Dept,  of  Mech.  Engrg. 

ATTN:  (Library  (Route  io 

Or.  P.  Wegener) 

New  Haven  10}  Conn.  (1) 

INDUSTRIAL  ORGAN  UA  IION^ 

Allied  Research  Associates 
ATTNi  Library  (Route  to 
Dr.  T.R.  Goodman) 

43  Leon  Street 

Boston  DfKass.  (1) 

Bell  Airosystem 
ATTNt  Library 
P.  0.  Box  1 

Buffalo  S,  N.  Yo  (  1) 


Boeing  Scientific  Res.  Labs. 
ATTNf  Research  Library 
■^.0.  Box  3981 

Seattle  2‘iy  Washington  (1) 

Chance-Vouqh t  Aircraft}  Inc. 
ATTNi  Library 

Da  1 ias  }  Texas  ( l) 

Convai r 

Forth  Worth  Division 
ATTNi  Library 

Forth  Worth  1}  Texas  (1) 

Convair  -  San  Diego 

ATTNt  Engineering  Library 

San  Diego  12}  Calif.  (1) 

Convair  Scientific  Res.  Lab. 
ATTNt  Library  (Route  to 
Chief }  Appl.  Res. ) 

P.  0.  Box  950 

San  Diegu  12,  Calif.  (l) 

Cornell  Aero.  Labs.,  Inc. 

ATTK't  Library 

Buffalo  21,  N^  Y.  (1) 

Douglas  Aircraft  Co. ,  Inc. 
ATTNt  Library 
027  Lapham  Street 
El  Segundo,  California  (1) 

Douglas  Aircraft  Co.,  Inc. 

ATTNt  Library 

3000  Ocean  Park  Blvd. 

Santa  Monica,  California  (l) 

Flight  Sciences  Laboratory 
ATTNt  Library  (Route  to  Dr. 

J.  Itenberg) 

1965  Sheridan  Avenue 
Buffalo  23,  N.  Y.  (1) 

General  Electric  Company 
Aerosciences  Laboratory  -  MSVD 
(Libiary) 

3750  "D"  Street 
Philadelphia  24,  Pa. 


(1) 
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General  Eltctrle  Company 
Resaarch  Labora'itory 
P.  0.  Box  1088 

Schanactady  S,  N.  Y.  (1) 

Grumman  Aircraft  Engr.  Corp* 
ATTNt  Library 

Bathpaga  L.  I.,  N.  Y.  (1) 

Hughat  Aircraft  Company 
Raaaarch  A  Oavalop.  Labs. 

ATTHf  Library 

Culver  City,  California  (1) 

Lockheed  Aircraft  Corp. 

ATTNi  Library 
P.  0.  Box  5B1 

Burbank,  California  (1) 

Lockheed  Aircraft 

Hlssila  Systems  Division 
ATTMi  Library 

Palo  Alto,  California  (1) 

The  Martin  Company 
ATTNi  Library 

Baltimore  3,  Maryland  (1} 

The  Martin  Company 
ATTNi  Library 

Denver,  Colorado  (1) 

McDonnell  Aircraft  Corp. 

ATTNt  Library 
P.  0.  Box  516 

St.  Louis  66,  Missouri  (1) 

Nexth  American  Aviation,  Inc. 
Missile  Division 
ATTNi  Library 
12214  Lakewood  Blvd. 

Downey,  California  (1) 

Northrop  Aircraft,  Inc. 

ATTNi  Library 

Hawthorne,  California  (1) 


Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California  (1) 

Republic  Aviation  Corpoiatior, 
ATTNt  Library 

Parmingdale,  L.  1.,  N.  Y.  (1) 

Uniteo  Aircraft  Corporation 
Research  Dept.  (Library 
400  Main  Screet 

East  Hax'tford  B,  Conn.  (1) 

Unified  Science  *£  ncletos,  Inc. 
ATTNi  f,  Nrluitc,',  Predldont 
926  S.  kT^oy*‘  ^arpway 
Paeeden*.  Callfornta  (1) 


